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Zero- to ultra-low-field nuclear magnetic resonance (ZULF NMR) provides a new regime for the measure¬ 
ment of nuclear spin-spin interactions free from effects of large magnetic fields, such as truncation of terms that 
do not commute with the Zeeman Hamiltonian. One such interaction, the magnetic dipole-dipole coupling, is 
a valuable source of spatial information in NMR, though many terms are unobservable in high-field NMR, and 
the coupling averages to zero under isotropic molecular tumbling. Under partial alignment, this information 
is retained in the form of so-called residual dipolar couplings. We report zero- to ultra-low-field NMR mea¬ 
surements of residual dipolar couplings in acetonitrile-2- 13 C aligned in stretched polyvinyl acetate gels. This 
represents the first investigation of dipolar couplings as a perturbation on the indirect spin-spin /-coupling in 
the absence of an applied magnetic field. As a consequence of working at zero magnetic field, we observe terms 
of the dipole-dipole coupling Hamiltonian that are invisible in conventional high-field NMR. This technique ex¬ 
pands the capabilities of zero- to ultra-low-field NMR and has potential applications in precision measurement 
of subtle physical interactions, chemical analysis, and characterization of local mesoscale structure in materials. 

PACS numbers: 82.56.-b, 82.56.Dj, 76.60.-k, 76.60.Jx 


Nuclear spin interactions are of substantial importance for 
many fields, including chemistry, quantum information pro¬ 
cessing, and precision measurement of fundamental symme¬ 
tries. The most common technique for measuring such in¬ 
teractions is nuclear magnetic resonance (NMR), typically in 
large magnetic fields in order to maximize signal via higher 
nuclear spin polarization and sensitivity of inductive detection 
m. However, the only terms of the spin-coupling Hamiltoni¬ 
ans that may be observed in high-field NMR are those that 
commute with the Zeeman Hamiltonian, which effectively 
truncates many interaction Hamiltonians that possess differ¬ 
ent symmetry. Recently, however, NMR experiments have 
been carried out in the opposite regime of very small magnetic 
fields GHl, taking advantage of advances in hyperpolariza¬ 
tion BID and new detection modalities EME), which of¬ 
fer a significant time savings compared to earlier field-cycling 
techniques fT7l fTSIl . In zero- to ultra-low-field NMR (ZULF 
NMR), the strongest interactions are the local spin-spin cou¬ 
plings, which involve coupling tensors that are of different 
symmetry from the Zeeman Hamiltonian and are many orders 
of magnitude smaller in amplitude, thus permitting the direct 
observation of nuclear spin interactions that vanish at high 
magnetic fields. Such terms that are only directly observable 
in the absence of large magnetic fields include the antisym¬ 
metric /-coupling (of importance for measurements of chiral¬ 
ity and parity violation), several terms of the direct dipole- 
dipole coupling, and a number of as-yet-unobserved exotic 
interactions such as those mediated by pseudoscalar (axion- 
like) bosons 021121. which would lead to anomalous spin- 
spin tensor couplings mm, most of which do not commute 
with the Zeeman Hamiltonian. 


As a proof of concept, we present in this Letter direct ob¬ 
servation of the untruncated residual dipolar coupling between 
nuclear spins in a weakly aligning environment in the absence 
of an external magnetic field. Dipolar couplings have long 
been used in high-field NMR to provide structural informa¬ 
tion in addition to the chemical shift. Previous work demon¬ 
strated zero-field /-spectroscopy of several systems for chem¬ 
ical analysis |[234l25ll . Additional information may also be ob¬ 
tained from zero-field NMR spectra via application of weak 
magnetic fields m. In the regime where dipole-dipole interac¬ 
tions can be treated as a perturbation to the /-coupling, zero- 
to ultra-low-field (ZULF) NMR allows sensitive measurement 
of the dipole-dipole coupling tensor. However, direct dipole- 
dipole couplings observed in solids are typically on the order 
of tens of kHz, substantially larger than /-couplings, and co¬ 
herence and population lifetimes are often too short for current 
ZULF methodology. Furthermore, all dipole-dipole coupling 
terms average to zero in isotropic liquids El- Smaller, scaled 
couplings are obtained by weakly aligning the molecule of in¬ 
terest in anisotropic media, such as liquid crystals 1261428) or 
stretched gels ll29l (30) . where molecular motion is partially 
restricted, yielding residual couplings. Such techniques have 
found widespread use in high-field NMR for structural mea¬ 
surements of proteins and small organic molecules MM- 
Here, we investigate the effects of residual dipole-dipole cou¬ 
plings (RDCs) on the zero-field spectrum of a model XA3 spin 
system: acetonitrile- 2 - I3 C (' ’CHiCN, where we have found 
it valid to neglect the 14 N spin due to self-decoupling aris¬ 
ing from fast quadrupolar relaxation [1351 ) aligned in stretched 
crosslinked polyvinyl acetate (PVAc) gels. 

Polyvinyl acetate (PVAc) polymer sticks containing be- 


2 



FIG. 1. a) Schematic illustration of the change in symmetry that oc¬ 
curs during uniaxial stretching of the gel environment due to swelling 
with acetonitrile-2- 13 C. The change in the order parameter is illus¬ 
trated by the different three-dimensional shapes, b) Energy level 
structure of a partially ordered 13 CH 3 spin system. F is the total 
spin angular momentum, K is the total proton angular momentum, 
and / CH is the one-bond 13 C - 'H 7-coupling. Solid arrows indicate 
allowed transitions, dashed lines indicate forbidden transitions. Note 
that for this system, both couplings are negative. 


tween 1-6% v/v divinyl adipate (DVA) crosslinker were pre¬ 
pared in 5 mm NMR tubes as described in the Supporting 
Information. Anisotropic gels were prepared by adding ace¬ 
tonitrile to the tubes and allowing the polymers to swell for 
2 weeks. Because the polymers were confined to the NMR 
tubes in which they were cast, swelling was uniaxial and 
equivalent to stretching along the axis of the NMR tube. A 
schematic representation of the process is shown in Fig. [I]' a). 

In order to maximize the ZULF NMR signal, the samples 
were prepared using labeled acetonitrile-2- 13 C to which was 
added 5% v/v deuterated acetonitrile for the purpose of high- 
held NMR characterization. The molecular order parameter 
ED EH for acetonitrile in the stretched gel environment was 
determined by analyzing the quadrupolar splitting of the deu¬ 
terium resonance ll26l[38 il using a 14.1 T NMR spectrometer 
with deuterium frequency 92.1 MHz (for additional details, 
see Supporting Information). The value for the electric-held 
gradient around the deuterium nuclei in acetonitrile was ob¬ 
tained from the literature f38l . 

The ZULF NMR apparatus has been described previously 
0151 [39|. Samples were pre-polarized in a 2 T permanent 
magnet located outside of the magnetic shielding for ~20 s 
and then shuttled pneumatically to the zero-held region over 
0.5 - Is. NMR signals were detected with an atomic magne¬ 
tometer featuring a 0.6 x 0.6 x 1.0 cm 3 87 Rb vapor cell op¬ 
erating at 180°C. Transient signals were collected over ~20 s. 
The spectra in Fig.[2]are the average of between 256 and 1024 
transients, and the spectra in Fig. [3] are the average of 8 tran¬ 
sients. 

The spin Hamiltonian in the presence of /-couplings and 


dipole-dipole interactions is 

y-{ = ft ^ J jk lj ■ L 

j;k> j 

-'^ZyrN 1 ' 0) 

j*>j Jk 

where ti is the reduced Planck constant, //q is the vacuum per¬ 
meability, jj and Yk are the gyromagnetic ratios of spins I,- 
and Ij., and r ,-j. is the internuclear vector connecting the spins. 
In the case of isotropic liquids, the dipole-dipole interaction 
term averages to zero. However, in aligned samples, such 
as stretched gels, the motional averaging of the dipole-dipole 
term is incomplete. For the system studied here, the swelling 
of the polymer gel with acetonitrile along the axis of the NMR 
tube leads to an orientational probability distribution of the 
solvent molecules that is slightly anisotropic, with the prefer¬ 
ential alignment axis (the director) determined by the swelling 
direction 12911301 . This axis is collinear with the sensitive di¬ 
rection of the detector, and is denoted z■ Because of the rapid 
rotation of the acetonitrile methyl group and the axial sym¬ 
metry of the alignment medium, the x and y components of 
the r j k vectors are averaged to zero, and the z components are 
scaled by the degree of alignment. Considering these averag¬ 
ing effects on the second term of Eq. <G the residual dipolar 
coupling Hamiltonian is 

'Hrdc - ~h ^ Djk (3/ jJk,- - 1/' U), (2) 

i\k>j 
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and where 9jk is the angle between the internuclear vector and 
the C 3 -axis of acetonitrile. We may also define the coupling as 
being directly proportional to the molecular order parameter 
•Szz = 2 (3 cos 2 Q z - 1 ) EH33, where 0- is the angle between 
the C 3 -axis of acetonitrile and the laboratory ’-axis and the 
time average provides a measure of the extent of alignment 
between the two axes systems, such that 


Djk - 


Fo yjVkli 1 


47T 


jk 


- (3 cos 2 (pjk - 1)5;.,, 


(4) 


where <pjk is the angle between and the C 3 axis, n/2 for 
Z>hh and the tetrahedral angle (2 arctan V2) for Dch- Because 
the order parameter S zz may be readily extracted from high- 
held deuterium NMR spectra, the only remaining parameters 
for the dipole-dipole coupling strengths are the angles and dis¬ 
tances that define the molecular geometry, which may be ob¬ 
tained from the literature l40l . 

It is worth pointing out the difference between the effective 
residual dipole coupling Hamiltonian in zero held and in high 
held. The zero-held heteronuclear coupling term in terms of 
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FIG. 2. Zero-field spectra of acetonitrile-2- 13 C with different degrees of ordering arising from the concentration of the cross-linker divinyl 
adipate (DVA). (a) K = I and (b) K = | peaks, (c) Peak positions as a function of molecular order parameter. The lines are calculated 
transition frequencies from Eqs. ( |8p4) with no free parameters. The order parameter for each sample was calculated from high-field deuterium 
quadrupole splittings using literature values for the electric field gradient >38) and the dipole-dipole coupling strengths were calculated from 
Eq. 0 using bond lengths from Ref. (4Ql . Solid symbols represent allowed A m F = 0 transitions, open symbols represent transitions with 
A m F = ± 1. Dashed lines indicate possible transitions that are not resolved. 


the total 'H angular momentum K and the l3 C angular mo¬ 
mentum S is 

rj-f(het) _ _^£) C h (3K-S - - K • S) , (5) 

as compared to the high-field case, 

, H< tru "c ) = -2HDchK z S : , (6) 

wherein a term of the form 

'Hrdc ~ ^H ( ' runc) = \d C k ( K + S _ + K.S + ), (7) 

is truncated because it does not commute with the high-field 
Zeeman Hamiltonian. The zero-field Hamiltonian is untrun¬ 
cated, and thus includes this so-called “heteronuclear flip- 
flop” term that is invisible to high-field NMR. 

In the regime where D,k <sc Jjk, the residual dipolar cou¬ 
plings may be treated as a perturbation on the /-coupling, 
yielding energy shifts that are calculated in the Supporting In¬ 
formation, and summarized schematically in Fig.[T]b), along 
with the allowed transitions. 

Because the observable in this experiment is the z- 
magnetization M z (t) oc Tr {p(t) YujljzYj)^ the detectable coher¬ 
ences are those with A F = 0, ±1 and A m F = 0. An addi¬ 
tional selection rule, AK = 0, arises in the case of equiva¬ 
lent spins (e.g. the methyl protons in acetonitrile) because K 2 
commutes with the Hamiltonian If4l l39ll . It follows that there 
is one allowed transition between K — \ states, between the 
|F = 0 ,m F = 0) and |F = 1 ,m F = 0) states. Based on the 
first-order energy shifts presented in the Supporting Informa¬ 
tion, this transition has frequency 

w o’o = ^ CH + Dch- (8) 


In addition, there are nominally forbidden (assuming that the 
detector is only sensitive in the ’-direction) A nip = ±1 transi¬ 
tions with frequency 


For the transitions between K — 3 states, there are allowed 
transitions with frequencies 

OJ 2 {° 0 — 2/ch + 2 D C h, (10) 

w t±i = 2 /ch + ^ (Dch + 3 D H h) , (11) 

and nominally forbidden transitions A trip = ±1 with frequen¬ 
cies 

w u±i = 2/ch - - (7 D C h + 3 D hh ) , (12) 

^ro* = 27 ch + - (5 D C h - 3 D HH ), (13) 

w U±i = 2 /ch + - (5 D C h + 9 D H h) ■ (14) 

If the detector axis is not exactly aligned with the direc¬ 

tor/quantization axis, the nominally forbidden transitions be¬ 
come observable. 

Zero-field spectra of acetonitrile-2- 13 C ( 13 CH 3 CN) in 
stretched polyvinyl acetate gels are shown in Fig. [2} a) for 
increasing values of the molecular order parameter. As the 
order parameter increases, the K = \ peaks corresponding to 
the ordered portion of the sample split, while the K — \ peak 
corresponding to excess isotropic liquid external to the gel re¬ 
mains unchanged. The lower-frequency peak in Fig. 0a) cor¬ 
responds to the A nif = 0 transition described by Eq. |8]( and 
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the higher-frequency peak corresponds to the A nip — +1 tran¬ 
sition described by Eq. The magnitude and phase of the 
A nip = ±1 peak are determined by the projection of the ini¬ 
tial spin-state population onto the transverse component of the 
detection operator, and is thus a signature of imperfections in 
the experimental configuration. Because the A mp = ±1 peaks 
are consistently narrower than the A mp = 0 peaks in Fig. [2ja), 
it appears that the linewidth is affected by inhomogeneity in 
the gel producing a distribution of order parameters and thus 
a spread in transition frequencies proportional to D (ii . 

Figure |2jb) shows four K = f peaks, three from the aligned 
acetonitrile-2- 13 C, and one from the isotropic liquid. The two 


lower-frequency peaks arise from the A mp = 0 transitions de¬ 
scribed by Eqs. ( flOpT) and the small higher-frequency peak 
corresponds to the A m F = ±1 transition described by Eq. 
Transitions corresponding to Eqs. ( T3p4| i are not resolved. 

In high-field NMR, terms in the Hamiltonian that do not 
commute with the Zeeman Hamiltonian are neglected, due to 
their immeasurably small effect on the NMR spectrum. Due 
to the absence of a large Zeeman interaction, ZULF NMR 
provides spectroscopic access to all spin-coupling terms HD . 
In traditional high-held NMR, only part of the heteronuclear 
dipolar coupling, '7Y (rn " !c) , f rom Eq. 0, yields measurable ef¬ 
fects in the spectmm. By itself, this term would yield no shift 
of the | F = 0 ,iiif = 0) «-» | F = 1 ,mp = 0) transition. In 
Fig. 0a), the residual dipolar couplings clearly shift the peak 
relative to the isotropic liquid showing the observation of the 
untruncated Hamiltonian of Eq. ([2]), including the contribution 
of Eq. 0, an interaction “invisible” to traditional NMR. The 
absence of truncation that permits observation of this term has 
also been demonstrated via the preparation of heteronuclear 
spin-singlet states in Ref. Il42l . 

As shown in Fig. 0c), the frequency shift varies linearly 
with the order parameter. The data closely match the simu¬ 
lated curves, which are calculated from the order parameter 
using Eqs. 0 and 0[T0. 

In the situation where the sensitive axis of the magnetome¬ 
ter is parallel to the director of the gel alignment, only the 
A nip — 0 transitions are detected. If the detection axis (or di¬ 
rection of alignment) is rotated, the A mp = ± 1 transitions are 
also observable, leading to peaks at higher frequency, with in¬ 
tensity dependent on the initial spin state populations and the 
angle between the detection and alignment axes. The addi¬ 
tional peaks appear in Fig.0because the measurements were 
carried out using a magnetometer configuration that featured a 
rotated axis of detection due to a non-zero effective field at the 
Rb cell. We attribute this effect to imperfections in the magne¬ 
tometer configuration, potentially related to laser alignment, 
AC Stark shifts, or a combination of the two. We point out 
that the effect was diminished after expanding the pump beam 
(thus decreasing the laser power density) and subsequently re¬ 
aligning the optics for the experiments in Fig. [3] It is worth 
noting that the rotation of the detection axis is necessary for 
the detection of the higher-frequency peaks only due to the 
axial symmetry of the system under study, which causes the 
terms of the dipolar coupling Hamiltonian that depend on the 
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FIG. 3. Acetonitrile-2- 13 C K = I peaks as a function of magnetic 
field applied orthogonal to the director axis, collinear with the mag¬ 
netometer pump beam. 


azimuthal angle to be averaged to zero. If these terms are not 
averaged to zero (e.g. in a biaxial phase), they will further lift 
degeneracy and yield additional peaks l43l . 

We have also investigated the effect of applied magnetic 
fields on the spectrum, as shown in Fig. [3] When the effective 
detection operator is collinear with the gel director axis, only 
the A nip = 0 transition is observed (corresponding here to a 
10 nT applied field). As the field is increased, however, rota¬ 
tion of the nuclear spins mixes states with different mp, and 
rotation of the alkali electron spins elicits a change in the sen¬ 
sitive axis of the detector. The overall result is that as the field 
is increased, the vectors defining the detection operator and 
the quantization axis cease to be collinear. This in turn leads 
to the observation of A nip = ±1 transitions, which become 
dominant above 50 nT, at which point the effective detection 
operator has been rotated substantially away from the director 
axis. 

In the regime where the Zeeman interaction strength is on 
the order of the residual dipolar coupling, the peak frequencies 
in Fig. [3] do not vary linearly with the applied field strength. 
This is because the dipole-dipole coupling Hamiltonian does 
not commute with the Zeeman Hamiltonian, and thus first- 
order perturbation is no longer sufficient to describe the sys¬ 
tem. 

In conclusion, we have demonstrated the direct influence 
of the heteronuclear dipole-dipole coupling “flip-flop” term 
(which is not directly observable in the high-held regime) on 
ZULF NMR spectra. The results are in agreement with a zero- 
free-parameter model of residual dipolar coupling utilizing lit¬ 
erature values for internuclear distances and the electric held 
gradient of acetonitrile, the latter being used to determine the 
order parameter via high-held deuterium NMR. The sub-Hz 
resolution of ZULF NMR may be of use for chemical and 
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structural analysis of small molecules, as well as for preci¬ 
sion measurement searches for anomalous spin-spin couplings 
along the lines of Ref. ED. 

In principle, all terms of the spin coupling Hamiltonian 
are observable in ZULF NMR, increasing the information 
available in NMR spectra. With appropriate systems (e.g. 
aligned chiral molecules) and continuing improvements in po¬ 
larization and magnetometer sensitivity, ZULF NMR may be 
a promising method to measure subtle interactions such as 
the as-yet-unobserved antisymmetric components of the J- 
coupling tensor ||44| . Based on calculations in Ref. |45l , mea¬ 
surement of a non-zero antisymmetric 7-coupling could per¬ 
mit the observation of a first-order energy shift arising from 
parity non-conservation in the molecular Hamiltonian. 

This technique may also find application as a probe of ma¬ 
terial structure, allowing for the measurement of interactions 
with lower than azimuthal symmetry. For example, measure¬ 
ment of the full anisotropic spin-spin coupling tensor of small 
molecular probes within porous materials may provide a more 
complete description of local geometry than do measurements 
of total surface area or average pore diameter. 
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